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I. 


INTRODUCTION 


The  continuing  search  for  lighter,  stronger,  more  economical 
missile  structures  has  led  to  the  investigation  of  fiber  reinforced 
composites  as  a possible  applicable  type  of  construction.  The  inherent 
properties  of  this  material  pose  many  difficult  problems  for  the  designer 
and  analyst.  Tests  on  actual  structures  have  proven  :hat  this  material 
is  unlike  conventional  structural  materials  in  many  respects.  Although 
much  has  been  learned  about  the  material  behavior  in  the  past  few  years, 
there  are  still  many  areas  that  are  unknown  and  unpredictable. 

An  initial  effort  [1]  used  a straight  line  element  for  a shell  of 
revolution  and  included  shear  deformations.  This  gave  very  good  results 
for  single-layered  shells. 

Nickell  and  Sato  [2]  used  a curved  shell  element  but  did  not  con- 
sider transverse  shear  deformations.  A single  layered  shell,  using  a 
curved  shell  element  and  considering  transverse  shear,  has  been  com- 
pleted [3].  Several  selected  numerical  examples  were  solved  and  the 
solutions  are  in  good  agreement  with  known  results. 

In  this  effort,  the  finite  element  method  is  used  with  a curved 
shell  element  considering  a laminated  shell  of  revolution  and  transverse 
shear  deformations.  The  field  equations  similar  to  Reissner's  theory 
of  thick  plates  [4]  are  used  as  a guideline  for  formulating  the  shear 
deformation  degree  of  freedom.  The  procedure  employed  is  similar  to 
that  of  Clough  and  Felippa  [5].  The  classical  Kirchof f-Love  assumption 
for  normals  to  the  midsurface  is  relaxed  in  favor  of  the  assumed  shear 
deformation  mode. 

II.  ELEMENT  GEOMETRY 

The  shell  to  be  considered  is  axisymmetric;  therefore,  it  is 
sufficient  to  define  only  the  shape  of  its  meridional  curve.  The  finite 
element  method  will  be  used  for  this  analysis.  The  element  is  shown  in 
Figure  1. 

For  this  analysis,  x - y are  the  local  rectilinear  coordinates, 
and  r - z are  the  global  coordinates.  The  angles  shown  in  Figure  1 
are  related  by  the  relation 

• + * + P - J (l) 
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Figure  1.  Curved  shell  element. 


The  angle  <t>  Is  the  angle  between  the  normal  to  the  reference 
surface  and  the  axis  of  revolution.  The  angle  ijr  is  the  angle  between 
the  chord  of  the  element  and  the  z-axis.  The  angle  p is  defined  to 
be  the  angle  between  the  chord  line  (the  x-axis)  and  the  tangent  to 
the  curved  surface. 


From  Equation  (1), 

sin  P » cos  (4>  + \|r)  = cos  $ cos  - sin  0 sin  ^ 
cos  P = sin  ($  + \|r)  = sin  <t>  cos  \|<  + cos  $ sin  \(( 


(2) 


To  approximate  the  meridional  curve,  the  following  substitute 
curve  is  assumed: 


(3) 


where  t ■ chord  length  of  the  element.  (0  < x < l) 


Differentiating  Equation  (3)  with  respect  to  x 

dy  2(-a2  ' al^  3^a3  ' a2^  2 4^34  ’ ®3^  3 5a4  4 

T"  = a.  + x + X x + i 7~  x 

dx  1 1 l2  l*  i4 


d2v  2^a2  ' al^  6^a3  " a2^  12^a4  " a3^  2 20a4  3 

72 — 7 — + — ^ — x + — ~s — x -~prx 

dx  l i l 


The  constants  a^,  a 2>  a^,  and  a^  can  be  determined  by  evaluating 
Equations  (4)  at  the  end  points 

a^  = tan 


„ Q l 3Q 

a2  = tan  P1  " 2R^  S6C  P1 


a3  = sec  P2  + R-  sec  ‘ 4 tan  ^2  " 5 tan  ^1 

a^  = ~~  sec362  - ^ sec3g1  + 3(tan  ^ + tan  P2>  . 

The  following  geometrical  relations  are  given  with  respect  to  the 
element: 

^r  r2  " rl 

Az  = 2^  - Z2 


^(Arl2  + (Az)2 


. Ar 

sin  V = — 


i Az 

cos  l * — 

After  the  substitute  curve  nas  been  established,  all  the  geometric 
quantities  can  be  written  as  follows: 


tan  P 


r » r^  + x sin  V + y cos  \)r 


(7) 


dr 

■^  = sin  \|r  + tan  P cos 


S " C°S  P £ 


d6  dBdS  1 

to  = dS  to  = ~ R S6C  P 


Since 


= 

dx 


tan  P 


fl  = to  (tan  P)  = Sec2p  to  ’ 

dx 


therefore 


dg  = _ 1 


sec  g 


dx  R 

The  quantity  dp/dS  is  negative  since  P is  decreasing  as  S is  increasing. 
Therefore 


Cx  = _ I 

a 2 R 

dx 


sec  p 


(8) 


and 


cos  0 = sin  P cos  i|r  + cos  P sin  ^ 
sin  $ = cos  P cos  i|r  - sin  P sin  \|r 

III.  TRANSFORMATION  OF  COORDINATE  SYSTEMS 


(9) 


The  displacement  vector  of  a material  point  on  the  midsurface 
in  the  local  principal  curvilinear  shell  coordinate  is  denoted  by 


{f  }T  - [u  , w , X , 7 ] 


(10) 
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where 


uc  = the  displacement  along  the  meridian. 


= the  transverse  (normal)  displacement. 


X = the  rotation  about  a meridional  tangent, 
c 


y = shear  deformation, 
c 


The  displacement  components  which  refer  to  the  local  rectilinear 
coordinates,  x-y,  are 


if  } » [U  , W , X , -v  ] 

r r’  r r r 

and  to  the  global  coordinates,  R-Z,  are 
(f}T  - [u,  w,  y] 

The  transformation  betx^een  these  components  can  be  seen  as 
follows : 

If  } = [q  ] (f  } 

1 * /■»  1 a- 


(ID 


(12) 


if  } = [q  ] X f } 
r ^r 


(13) 


where 


[qcl 


cos  p sin  f 0 0 

-sin?  cos  P 0 0 

0 0 10 

0 0 0 1 


and 


(14) 


[qrl 


sin  4'  -cos  400 
cos  4'  sin  4 0 0 

0 0 10 

0 0 0 1 


(15) 
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IV. 


STRAIN-DISPLACEMENT  RELATIONS 

The  strain-displacement  relations  can  be  written  as 


du  w 
c c 

p a — ■■■  4*  

1 dS  R 


^ (u£  cos  + w£  sin  <t>) 

_d_  K ^ \ 

1 * " dS  \dS  R y 

cos  $ / dwc  uc\ 

K2  r ^dS  ' R / 


7ly  “ “r 

The  strains  defined  in  Equations  (16)  are  now  transformed  into  the 


local  rectilinear  coordinates  as  follows  (recall  that  dS  = 

e-  -|)  = 

du  „ dw 

e,  * “t”  cos2p  + -r^-  cos  p sin  P 
1 dx  dx 


cos  p 


e_  * — (u  sin  ilr  + w cos  i(r) 

2 r r r 

.2 

3„  wr  2 cos  P sin 

*1  * -cos  P TT £ 

dx  1 


. 2e  2 du  dy 

. sin  P - cos  P r r . _ _ „ 

+ ■ ■ r-_ — — - — — COS  P 

R,  dx  dx 


dw  „ d2u 

+ sin  p cos  p ~ 

dx  .2 

dx 


* / o dw  du  \ 

js  G / 2 „ r , „ „ r\  cos  $ 

(cos  P — — - sin  p cos  P — — ) - v 

r y dx  dx  J r r 


7iy  - -rr 
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In  matrix  notation  this  can  be  written  as: 


0 0 0 


0 0 


0 1 x 


0 0 


2 2 2 2 

0 -sc  0 c 2xc  3x  c 0 0 


0 0 0 


1 x 


where 

s = sin  p,  c = cos  p 
This  can  be  written  symbolically  as 

{fr)  = [♦]  {a}  (2: 

where  {o}  is  the  generalized  coordinates  vector  for  the  curved  shell 
element. 

The  shell  displacements  shown  in  Equation  (20)  represent  4 degrees 
of  freedom  at  a node,  two  translation,  two  rotation.  The  8 degrees  of 
freedom  connected  with  the  nodes  of  the  element  are  written  as  the 
displacement  vector 

{6r)  - [url,  wrl,  Xrl,  7rl.  ur2’  wr2,  Xr2,  rr2]  . (2, 

The  generalized  displacements  {a}  are  related  to  the  nodal  point 
displacement  vector  I&6}  by 

(cj  = [A  1 (6e)  . (2! 

r r 

la}  is  evaluated  as  follows: 

1)  At  x = 0 


“L  = utl 


CL  * W 

u3  rl 

a_  - y , 

7 rl 


-0^  sin  P^  cos  P1  + cos  p^^  = 


Therefore,  from  (23) 


where 


tan  p 


1 


*1  “ l 

2 tan  p^  + tan  P2 


a2  = 


a3  = 


tan  P^  + tan  P2 


1 


bl  = 2 


cos  P, 


2 2 

l COS  p. 


J i coa  p 


r 


£ cos  p„ 


5 2 

l COS  P0 


The  transformation  of  ib  )6  to  the  global  coordinate  lb}6  is  given 


bv 


Substituting  (25)  into  (23)  gives 


[A  ] <6  }e  = [A  1 [R]  (&}e  = [A]  {&)' 


where 


[Ar]  [R] 


-cos  V 


-a^  sin 


cos 


a,  sin  V -a.  cos 


a2  sin  i 

-3  cos  , 
.2 


-a2  cos  V 

-3  sin  ■» 
.2 


-a2  sin  ♦ a2  cos  » 
+ 3 cos  y + 3 sin  v 


-a^  sin  . 

+ 2 cos  , 
,3 


•f  2 sin  . 
.3 


a^  sin  U 

- 2 cos 
,3 


-a^  cos  v 


VI.  STRESS-STRAIN  RELATIONS 

For  an  axisymmetric  shell  of  revolution  subjected  to  axisym- 
metric  loadings,  the  stress  resultants  and  couples  can  be  expressed  as 


”2 

Ql)  . 


11 

E12 

E13 

E14 

0 

21 

E22 

E23 

E24 

0 

31 

E32 

E33 

E34 

0 

41 

E42 

E43 

E44 

0 

0 

0 

0 

E 

L L 


The  quantities  are  related  to  the  principal  curvilinear  coordinate 
system  with  1 as  meridional  direction  and  2 as  circumferential  direction 
Symbolically  this  can  be  written  as 


The  detail  derivation  of  [E]  is  given 


where  [E]  is  the  elasticity  matrix 
in  Appendix  A. 


Substituting  (19)  into  (17)  gives 


]T  (qrl  >P>  =< 


2 2 
p sin  '-  + p cos  V sin  + p cos  w - P sin  cos  \|/ 
rx  ry  *x  *y 

2 2 
Px  x sin  ; + py  x cos  v sin  V + px  x cos  \|»  - Py  x cos  i sin  V 

p cos  v sin  » + p cos^  v - P cos  V sin  v + p sin^  v 

rx  y *x  ry 

2 2 
p x cos  *„  sin  + p x cos  \ - p x sin  V cos  v + p sin  * 
rx  ry  rx  ry 

2 2 2 2 2 2 
p x sin  ; cos  . + p x cos  - - p x sin  * cos  v + p x sin  U 

rx  y *x  ry 

3,  . . 32  3.  . ^ 3 2 

p x sin  . cos  - + p x cos  \ - p x sin  v cos  ^ + p x sin  „ 

x v rx  y 


)(45) 


r *\ 

p„ 


p x 
rx 


P x 

y 


p x 

y 


p x 

y 


(46) 


VII.  RESULTS  AND  CONCLUSIONS 


Various  appropriate  structures  were  analyzed  by  the  computer 
code  for  this  development.  Comparisons  were  made  with  known  solutions 
for  single  layered  structures.  The  results  were  extremely  close. 
Effort  is  now  being  made  to  find  some  results  for  laminated  structures 
so  that  a comparison  can  be  made. 


19 


The  accuracy  obtained  by  this  method  depends  directly  on  the 
extent  to  which  the  assumed  displacement  patterns  are  able  to  repro- 
duce the  deformation  actually  developed  within  the  element.  Since  the 
chosen  displacement  patterns  satisfy  the  requirements  of  completeness 
and  conformity  (continuity  of  displacement  at  element  boundary)  as  the 
size  of  the  element  decreases  indefinitely,  the  solution  obtained  from 
the  minimization  of  potential  energy  converges  to  the  exact  solution. 

There  still  remains  a need  to  add  geometric  and  material  nonlinear 
ties  to  this  analysis.  The  material  quickly  becomes  nonlinear  as  the 
matrix  material  begins  to  crack  or  "craze"  while  the  fibers  are  still 
intact.  Also  the  capability  for  buckling  prediction  is  one  of  the 
major  needs  for  this  type  construction.  These  items  are  now  being 
investigated  and  will  be  the  subject  of  a later  report. 


Appendix  A. 

ELASTICITY  MATRIX 


Individual  curved  finite  elements  can,  in  general,  be  composed  of 
a number  of  anisotropic  layers  of  varying  thickness  along  the  meridional 
coordinate.  For  a single  lamina,  considering  shear  deformations,  the 
constitutive  relation  is  given  as 


°L 

O' 

nl 

q12 

0 

0 

0 

€L 

aT 

O' 

V12 

Q22 

0 

0 

0 

€T 

tlt 

0 

0 

Q44 

0 

0 

7LT 

V 

0 

0 

0 

Q55 

0 

7h 

V 

0 

0 

0 

0 

V 

7T£ 

where  the  transverse  normal  stress  cr^  has  been  omitted  and  the  laminae 

are  orthotropic  with  respect  to  the  principal  elastic  axes  L-T.  These 
axes  need  not  coincide  with  the  axes  of  the  curvilinear  coordinate 
system  1-2,  (Figure  A-l),  (1  is  the  meridional  direction)  and 


O'  . 
W11 

V(1  - 

VLT  VTL* 

O'  . 
w12 

v E / 

lt  r 

(1  ' VLT  VTL 

m 

v E / 

'tl  1/ 

(1  ' VTL  VLT 

O'  - 
^22 

ET/d  - 

VLT  VTL) 

Q44- 

glt 

Q55- 

V 

Q66 

• 

(A -2) 
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Figure  A-l.  Material  axes. 

Equation  (A-l)  can  also  be  written  for  the  Kth  layer  in  the  following 
forms : 


(A-3) 


!f 


Qll  Qu  cos  3 + 2(Q,o  + 2 Q/./, ) sin  3 COS  9 + Q;o  sin  ) 


‘12  x44 


^12  = ^11  + ^22  ' 4 Q44)  sin2  3 cos2  3 + Q^2  (sin4  9 + cos4  9) 

Q22  = sin^  1 + 2(Q|^  + 2 ^44  ^ Sln2  4 cos2  9 + ^22  cos ^ 5 (A-5) 

= (Q[i  + Qu  - 2 ^4>  sin  9 cos3  9 + (Q[2  - Q22  + 2 Q44>  sin3  9 cos  9 


Q24  = (Qll  ' Q12  ‘ 2 Q44}  Sin  3 cos  9 + (Q|0  - Q;o  + 2 Q'A)  sin  9 cosJ  9 


12  x22  x44; 


Q44  “ + Q22  * 2 Q^2  ’ 2 ^44  ^ Sin2  9 cos2  9 + ^44  (sin4  3 cos4  9) 

*55  = Q55 
^66  = Q66 

Substituting  the  midsurface  strain  and  curvatures  into  Equations  (A-4) 
the  following  expression  is  obtained: 


ao  ) = IQ]  Je„  +{  [Q]  , < 


(A-6  ) 


By  integrating  over  the  total  thickness  of  the  laminate,  the 
generalized  stress  resultants  in  terms  of  midsurface  strain  and  curva- 
ture are  given  as 


[SI  = Z [QW1  (hk  - hk_j) 


in  which  hfc  and  h^  ^ = the  distances,  respectively,  from  the  midsurface 
to  the  inner  and  outer  surfaces  of  the  k-th  layer. 


Appendix  B. 

ELEMENT  STIFFNESS  MATRIX 


The  element  stiffness  matrix  is  found  by  writing  the  total  poten- 
tial energy  of  the  axisymmetric  shell  of  revolution  and  minimizing  it 
for  the  imposed  constraints  and  loading  conditions. 

The  potential  energy  for  a linear  elastic  shell  of  revolution  in 
the  absence  of  thermal  and  body  forces  can  be  formulated  as  follows: 

///  I {a}  dv  - I (f}T  {P}  dA  (B-l) 

V A! 

where  the  vectors  (e),  {a}.  If},  and  {p}  represent  the  strain,  stress, 
displacement,  and  equivalent  surface  traction  vectors,  respectively. 

Introducing  the  stress  resultant  vector 

{S>  = t {0}  (B-2 ) 

where  t is  the  thickness  of  the  shell.  Equation  (B-l)  may  be  written  as 

n = 17'  i (e}T  Is)  — - /;  (f}T  {p}  dA  . (B-3) 

v 2 c \ 

The  first  integral  is  evaluated  over  the  entire  volume  V of  the 
shell  and  the  second  over  the  portion  A^  of  the  midsurface  of  the  shell, 

where  the  equivalent  surface  tractions  are  prescribed.  Since  the  state 
of  displacement  throughout  the  shell  is  defined  element  by  element,  the 
total  potential  energy  may  be  considered  as  the  sum  of  the  potential 
energies  of  all  individual  elements,  i.e.. 


The  potential  energy  contribution  of  element  "e"  will  now  be 
considered.  The  state  of  displacement  defined  for  the  element  in 
local  rectilinear  coordinates  x-y  can  be  expressed  in  matrix  form 
in  Equation  (21)  as 

(fr)  = (01  {a}  = [0]{Ar]  {5®}  . (B-4) 


, * PRECISING  rJ£E  BLANK- NOT  FILMED 


Transformation  of  {fj  into  the  global  coordinate  system  may  be 
obtained  from  Equation  (13) 


if)  = [q  ]T  {f  } = [n]  {5  }< 
r r r 


(B-5: 


where 


[N]  = [qr]T  [*]  [A  j 


(b-6; 


“V  e COlU”  vect<>r  V tta  eight  discrete  parameters 

°f  the  '>“*«  *•  given  in  Equation  (25b) 
def ined^isplacement  ^pattern.  ^ L.riX 

diepi^^ItSIre  <«>  -rein. 


= [B]  (b}‘ 


where 


(B-7 


mi 


I'*"  J l A ] 


to“t£e“  rei1>vect«"atS  Ir'!^"8  !he  »°1"t  ‘■‘•H-c-.nt  vector 

expressed  as  8 “C  stress-^rain  relations  can  be 


<s)  = [E]  {«■} 

(B-9) 

where  [E]  is  a function  of  the  elastic  properties  of  th*  .in 
Each  element  can  be  assigned  different-  f th  element- 

relations  in  Equations  (B-91  Jr a S « properties.  If  the 
(B-3 ) the  norenft  1 ( h (B'5)  and  are  substituted  into 

1 he  P°tential  energv  contribution  for  the  element  becomes 

- in  i (e'J1  IB)1  |E)  |B,  lb'}  f . ff  hVT  WI  (p)  (B_ 


where  V,  la  the  <.lu.e  of  the  element  end  A,  is  that  pert  of  the  mid. 

surface  area  of  the  element  which  coincides  with  the  midsnrface  area  A 
of  th.  shell  over  which  the  equivalent  surface  tractions  are  prescribed 
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Since  the  discrete  parameters  fe®}  are  not  a function  of  spatial 
coordinates,  the  potential  energy  of  the  element  may  be  written  as 


it®  = {5e}T  ///  i [B]T  [E] [B]  f 


{&®}  - {&e}T  fl 
A1 


ff  [N]T  {p}  dA.  (B-ll) 


Since  the  assumed  displacement  patterns  for  each  element  satisfy 
various  requirements  such  as  completeness  and  conformity,  the  best 
values  that  can  be  obtained  for  the  total  nodal  point  displacements 
of  the  finite  element  representation  of  shells  of  revolution  are  those 
that  minimize  the  total  potential  energy  of  the  shell  under  the  con- 
straints imposed;  i.e. , the  best  value  of  15/  are  those  that  satisfy 
the  system  of  linear  equations 


(B-12) 


where  (6)  is  the  total  nodal  displacement  vector  of  the  system. 

In  forming  the  system  of  Equations  (B-12),  it  is  convenient  to 
have  an  expression  for  the  spatial  derivatives  of  the  potential  energy 
of  each  element  "e"  with  respect  to  its  own  nodal  point  displacement 

vector  {b6} , i.e., 


a {5e} 


T C'W®  on6  Sn®  5n®  c>«e  dn®  5n®  Sit® 
'WI  J’XI  \ ^ °WJ  5XJ 


(B- 13) 


By  use  of  Equation  (B-10),  this  expression  can  be  obtained  as 

—7 - rrr  [b]t  ie] [si  ^ {&e>  - 11  [n]t  {p}  dA 

5 {6e}  V C A, 

_ e j 1 1 

L e -J 


(B-14) 


The  terms  in  the  first  and  second  brackets  are  normally  defined  as 

0 

the  element  stiffness  matrix  [K  ] and  the  element  generalized  nodal 
point  force  {f6},  respectively.  Pence, 


[K6]  = ///  [B]T  [E] [B]  ^ 
V 

e 

{f®}  = //  [N]T  {p}  dA 


(B-15) 


(B-16) 
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By  properly  combining  the  submatrices  in  Equation  (B-14)  obtained 
for  each  element,  the  total  matrix  equation  representing  Equation  (B-12) 
can  be  constructed  as 

[K]  {b}  = {f}  (B-17) 

and  then  solved  for  the  nodal  point  displacements.  Once  the  nodal  point 
displacements  are  obtained,  the  corresponding  stress  resultants,  stresses, 
and  strains  for  the  defined  displacement  patterns  can  be  calculated  from 
Equations  (B-7)  and  (B-9). 

If  Equation  (B-8)  is  substituted  into  Equation  (B-15)  and  the 
volume  increment  for  a shell  of  revolution  is  taken  as 

dV  = 2«t  R(x)  dx,  (B-18) 

then  the  element  stiffness  matrix  for  the  axisymmetric  shell  element 
takes  the  form 

IK®]  = 2n  / [B]T  [El  [B1  R(x)  dx 
o 

= 2 x [A]T  [G]  [A]  (B-19) 

where 

[G]  - J[*']T  [E]  [0 ’ 1 R(x)  dx  . (B-20) 

o 

The  Integration  is  over  the  chord  length  of  the  meridian  cross 
section  of  element. 

It  is  assumed  that  the  equivalent  surface  traction  over  the  mid- 
surface area  A^  where  tractions  are  prescribed  varies  linearly  between 
the  two  nodal  circles  I and  J.  That  is, 

(p  }T  = [0  (P  + P’  x)  0 0]  (B-2 1) 

c n n 

where  (p  } is  the  surface  traction  vector  expressed  in  local  curvilinear 
c 

coordinates.  Transforming  into  global  coordinates  the  following  is 
obtained: 

{ P > “ [q_lT  (qJT  <P  > • (B-22 ) 

r c c 
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Substituting  Equations  (B-6)  and  (B-22)  into  Equation  (B-16)  the 
generalized  element  nodal  force  vector  becomes 

If6}  = 2 it  |VrIT  MT  [qrl  [qrlT  [qc]T  (pJ  R(x)  dx 

or 

(Fe>  = 2«[ArlT  f [*]T  [qclT  (Pc)  R(x)  dx 


-sin  p 
-x  sin  P 


cos  P 


‘c  c n 


x cos  P 


x cos  P 


-x  sin  P 
-x^  sin  P 


x cos  P 


T T r 1 X COS  p 1 I X COS  p 

M [qj  (p J = P V \+  P'  l 


X COS  P 
4 

X COS  P 


1. 
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